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For a generic f(R) which admits a polynomial expansion of at least third order (i.e. -j-J^ 7^ 0) 
we find the near-throat wormhole solution. Necessary conditions for the existence of wormholes in 
such f(R) theories are derived for both zero and non-zero matter sources. A particular choice of 
energy-momentum reveals that the wormhole geometry satisfies the weak energy condition (WEC). 
For a range of parameters even the strong energy condition (SEC) is shown to be satisfied. 



T—l 1 I. INTRODUCTION 

o: 

Black hole solutions in modified theories of gravity (namely f(R) gravity) has received attentions in recent years 
[![. Parallel to black holes there are studies on wormhole solutions in f(R) gravity as well 0j- Upon derivation of the 
necessary conditions for the existence of black holes in f(R) gravity [3[, it will be in order to explore similar conditions 
for wormholes in the same theory. This amounts to make series expansions for functions around the throat radius 
I 1 (r = r ) of the wormhole. That is, any function h(r) can be expressed as h(r) = h(r ) + h'(r )(r — r ) + 0((r — r ) 2 ). 
Such expansions are carried out for the Ricci scalar R, f(R), F = ^ and higher derivatives, as well as the pressure 
and density of the energy-momentum tensor. Let us remind that even in the absence of external sources the curvature 
sources of f(R) theory remain intact. The reason for resorting these expansions is entirely physical: near the throat 
O"!' radius of a wormhole designated by r — tq = x > 0, (x 2 < 1) an observer doesn't feel anything unusual so that the 
5^ \ functions specifying all salient features are expressible in analytic functions. Once these expansions are substituted into 
(DJQ Einstein's equations and zeroth, first and higher order terms in | r — ro | arc identified automatic constraints emerge. 

These will comprise the necessary conditions for the existence of wormhole throats in any viable f(R) theory as a 
^v^j . result of such a " near-throat test" . At the zeroth order the conditions are easily tractable and yield tangible results 
" to say enough about an f(R) whether it admits a throat or not. We must note that the higher order terms, even the 
I/"") , first one turns out to yield rather complicated relations in implementing the test. In this study we concentrate on an 

,3 r 

f(R) theory that admits at least a non-zero third order derivative, i.e. 7^ 0, while all the rest will be dictated by 
Einstein's equations. In this study our ansatz metric involves two unknown metric functions to be determined: the 
redshift function $ (r) and the shape function b (r). Upon expansions aforementioned in powers of (r — ro) there will 
be severe restrictions on these functions. The simpler case consists of choosing $ (r) as constants, however, herein this 
will not be our strategy. One thing observed is that in the absence of matter sources the weak energy condition (WEC) 
is violated in the construction of wormholes. To overcome this problem we introduce external energy-momentum and 
• • ■ search for the satisfaction of the energy conditions. We find that WEC holds true while the strong energy condition 
(SEC) becomes valid under more stringent regulations. Integration of f(R) incorporates constants that are related 
to the cosmological constant and with / without that constant validity of our energy conditions remains still valid. 
After we determine the redshift and shape functions analytically we embed our geometry into the (r, z, (f>) sector with 

the embedding function z = irocosh" 1 C^j • The shape of our traversable wormhole is depicted in Fig. I for the 
particular throat radius ro = 1. 

The paper is organized as follows. In Sec. II we introduce our formalism of expansion around the throat in a 
generic f(R) theory. Section III extends / generalizes the formalism in the presence of external matter sources. With 
Conclusion in Section IV we complete the paper. 
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II. THE FORMALISM 



We start with a general action for f(R) gravity written as 



S=±- I V=9f(R)d*x (1) 



2 k 



in which n = 8irG = 1, and / (R) is a real arbitrary function of the Ricci scalar R. Obviously, in the presence of 
physical sources the action is to be supplemented by a matter part Sm- The 4— dimensional standard form of the 
wormhole line element in spherical symmetry is given by Q 

where <£> (r) and b (r) are the redshift and shape functions, respectively. The throat of the wormhole is located at 
r = ro, at which b (ro) = ro and radial coordinate r is larger than ro. We note that $ (r) and b (r) , in addition to R, 
f(R) and its higher derivatives should also satisfy constraints to have a traversable wormhole. These conditions for 
b(r) and $ (r) are: i) b' (ro) < 1, ii) b (r) < r for r > ro and iii) e 2 *( r ) must not have any root (horizon) i.e., $ (r) 
must be finite everywhere. 

Variation of the action with respect to the metric yields the field equation 



FR^-t S »-W^F + 8»nF = (3) 



in which □ = V M V M = -^^(7=^) and VV^/i = g x "V\h^ = g x » (d x h lfi - T^hp) . The field equations 
explicitly read as 

FR\ - - + UF = V*V t F (4) 

FR r r -L+UF = V r V r F (5) 

FR e - { + UF = ^ 9 ^eF (6) 

which are independent. Note that the ip(p equation is identical with 89 equation. By adding the four equations (i.e., 
tt, rr, 99 and (ftp) we find 

FR-2/+ ZDF = (7) 

which is the trace of Eq. (3). 

Our method to solve these equations is as follows: First, we consider the Ricci scalar as a series about the throat 
i.e., 

R = Ro + R Q x + -R^x 2 + .... (8) 

in which x — r — ro > 0. Herein and in the sequel a prime stands for the derivative with respect to r and a sub 
implies that such a quantity is evaluated at the throat. Next, we expand all other functions involved in the field 
equations i.e., 

$ = $ + $ x + i$ V + ... (9) 
b = b + b' x + -b'^x 2 + ... (10) 
/ = /o + /oZ + i/oV + ... (11) 
F = ?L=F + F>x+ 1 -FZx 2 + .., (12) 

E = ^=E Q + E' x+ l -E''x 2 + ... (13) 

d 3 f 1 
H = dR^ = H ° + H '° X + 2 H ° x2 + - (14) 
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and finally by equating different powers of x we find all coefficients in terms of ro, Ro, R' , .... In zeroth order we find 

fa - r, K = i M* -i), *S = -j^jjy* - ^ <«> 

p 2r / _ R' Q r Q f _ r / 

^ " >gi? -2)i? ' /o -^R^2'^-" 2 ^^2' (16) 
while to first order in x we obtain 

°o - (#o r o + r o^o - 8J , $ — — -2 , (17) 

3r rg (r2i? - 4) 

u _ Jo(R'oro + R'o) F , _ /o(i%ro + i%) ,„ _ ^o/o (gcfo - 2fl ) 
iio - 2 <( r 2 il . _ 2) ^o-2 jR , 2(ro2i?o _ 2) , /o - ro2i?o _ 2 • 

We note that in our solution <f> an d /o are n °t specified. As one can see from the line element $q can be absorbed 
into the redefinition of time t, so we set it to be zero. On the other hand fo is a principal constant which we leave 
free so that other constants can be expressed in terms of it i.e., Q = r 'i^°_ 2 an d so on. Clearly this is just the zeroth 
order approximation and it can not be considered for the generic solution. 

In order to find the conditions b (r) must satisfy we start with b' (ro) < 1, which implies 

r 2 Q R < 4. (19) 
The other condition i.e., b (r) < r for r > ro is automatically satisfied since 

i( r gfl -l) + I^i + ... <l (20) 

and for small x it leads to rg-Ro < 4 which is assumed valid. The last condition which constrains $ (r) to be finite 
takes the form 



+ 2 _ rgijg + 8 - 2r 2 R + 2r 3 R' ^ 
ro(4~r 2 R ) X 2r 2 (r 2 R - if 



$( r ) = , 0JH) -ro !i0^0^ ^2 + ^ (21) 



This is guaranteed if x is small enough to keep ^ * 2R * finite together with rgi? < 4. 

Now let's consider the energy conditions. If we directly write the field equations in the form of standard Einstein 
equation 

G y » = ^TZ + fZ (22) 



in which Gf, stands for the Einstein's tensor, with 



rpV ^ 



V y V M F - OF - -/ + -i?F <5 



(23) 



and is the external stress energy tensor which in the present case is zero. Next, we consider the energy density p 
and the pressure components produced by the geometry of f(R) gravity in the form [B| 

f; = diag [-p (r) , Pr (r) , Pt (r) , p t (r)] . (24) 

After some manipulation we find at the throat (i.e. from the zeroth order Einstein's equations) 

P =-f° = r l^f±, (25) 



/■- 

and 



Vr=f\ = -\ (26) 



p e =p v =T 2 = —-^. (27) 



1-rjRo 

r 2 

It is observed that imposing 1 < r^Rg < 4, makes p positive but p§ — p v remain negative and in any case p r is negative. 
Note that it may be this negative pressure that protects the wormhole against collapse. From the summation of (25) 
and (26) it can be seen easily that the weak energy condition (WEC) which says that p+Pi ^ and p ^ 0, is violated. 
In the next section we may avert this situation by adding external sources. 



4 



A. Wormhole supported by a matter source 

In this section we consider a matter source which provides an energy-momentum tensor of the form 

= diag[-p,p,q,q] (28) 
in which p, p and q are arbitrary functions of r. The held equations, read now 

g^ = |t; + t;. (29) 

The same method as we used in the previous section, i.e., expansion of quantities about the throat, including p, p 
and q 

P = Po + Po x + - (30) 

P = Pa +Pa x + ••• (31) 

Q = qa + %x + ■■■ (32) 

would lead to the following results in the zeroth order: 

P °-2 fo+ ^^ 2^ ' (33) 



1 F (b{,-l)(^r + 2) 
Po--^fo + ^ (34) 



and 

l f , F % + 1) F '(6' - l) 
2 J ° + 2r§ 2r 



Qo- -77/0 + ^3 ^ • (35) 



Next, we try to apply WEC which implies p > and p +Pio > 0. One finds the following conditions which should 
be satisfied simultaneously: 

P >0^ h - F ^- b °K {1 - b '^° F °>0, (36) 
ro ro 

Po+Po>0^ r F^ + 2F < 0. (37) 

and 

l-b' 



Po + qo > -> F 



1 + r $ 



l + & 



> (38) 



in which we used the choice that b' Q — 1 < 0. Next, one may add to the foregoing conditions p + po + 2go > which 
implies the strong energy conditions (SEC) and gives 

2Sf£ + S£3>_ A > a (39) 

In the case of = [j| one obtains 

WEC -> f > — — ; — < ; > (and consequently Fq < 0) ; (40) 

F) nj 



and 

^0 F>r " r§ 2F r 



SFC -> WFC and < ^ - v °> < =2 - ov " o; . (41) 



We note that although these conditions are very complicated in principle one can choose the proper parameters to 
satisfy them. In the next section we give an exact solution to the problem and with a particular example we shall 
justify our prediction. 
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III. GENERALIZATION 



In this section we study the field equations not at the throat but at any r ro . This is also a generalization of Ref . 
[j| in which $ =constant. The field equations are given in Eq.s (4)- (7) and the line element is considered as (2) while 
the energy momentum is given by (28). A detailed calculation would imply the following relations for the components 
of energy momentum: 



P = 



f | 2r (*" + $' 2 ) (r - b) + [-36 + (4 - b>) r] <&' p + (6' - 4) r + 36 ^,, _ ^ _ b^ fl , 



2r 2 



2r 2 



(42) 



/ 2r 2 (b - r) ($' 2 + $") + (r$' + 2) (6V - 6) ( r - 6) (r$' + 2) , 
P = ~ 2 + 2^ F + ^ F ' 



/ 2r (6 - r) $' + r6' + 6 [-2r (r - b) + (&' - 2) r + 6] , 
Q ~ ~2 + 2^3 F 2T 2 

The case of Ref. [f| is easily observed if one sets $ =constant and therefore 



f , (y-4)r + 3b F , 
^=2+ 275 



1 - - | F" 

r . 



1 - - ] F". 



(43) 
(44) 

(45) 



/ (b'r-b) 2(r-b) . 



(46) 



/ ^ + (6 , -2)r + b F , | / ft 



2r 3 



2r 2 



(47) 



The other limit can be checked for r = r (b = bo) and the results found in (33)-(35) are recovered. 

Next, let's consider an isotropic velocity distribution which implies p = q together with the ansaetze b(r) = r^/r 

and $ (r) = 2 In ( — J . These in turn lead to 



r 2 (r 2 - r 2 ) F" + r (2r 2 - r 2 ) F' + 4r 2 F = 0. 



(48) 



This equation admits two independent solutions and to keep our calculation analytic, we choose the simpler one given 
I'.v 



in which C\ is an integration constant. Same ansatz in the Ricci scalar gives 

6 (rl - 2r 2 ) 



R = 



(49) 



(50) 



which implies (through F = -J^ = ^j) 



f 



24Ci 



2 \ 5/2 



1-^ 



C 2 



(51) 



for the integration constant C2 that is related with the cosmological constant. Finally one obtains 

\ 5/2 



/OR) 



24Ci 



( 



\ 



6 



c 2 



(52) 
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We note that the Ricci scalar R satisfies — §■ < R < 0. This causes no restriction on the above expression and is well 
defined everywhere for r > tq. Our next step is to find the stress-energy components which are as follow: 



(24r 6 - 12r 4 r 2 - 18r 2 r 4 + 6rg) d + 5r 5 r 2 C 2 ^r 2 - r' 2 . 
p = m 2 I r _ T == f ' V 53 ) 

(24r 6 - 52r 4 r 2 + 32r 2 r 4 - 4r%) C x + hr>r 2 C 2 ^/r 2 - r 2 
10 r o r V r — r o 

Furthermore the WEC implies that p > and p + p > 0, which in closed form read 

p > (24r 6 - 12r 4 r 2 - 18r 2 r 4 + 6r^) Ci + br b rlC 2 \Jr 2 - r 2 > (55) 

and 

d (4r 4 - 5r 2 r 2 + r 4 ) , 
r°y/r z — Tq 

One can see that both conditions are satisfied if both C\ and C 2 remain positive. In addition to WEC it is remarkable 
that the strong energy condition (SEC) i.e. p + p > and p + 3p > 0, is also satisfied under the following condition: 



C 2 3 3 - 16£ 2 + 8£ 4 



y/e^l. (57) 



Ci " 5 r 2 C 5 

Here C\ and C2 both are positive as chosen above for WEC and £ ( = — I is a positive parameter such that 1 < £ < 



v / 4+v / io ^ j t gjjQ^j |-, e added that even for the case of C 2 = 0, the SEC holds true in the particular range of £. Let 
add that an expansion of f(R) in powers of R yields 



OA.C 

f(R)°i '+C 2 + C 1 R + 0(R 2 ) (58) 
or 



which gives the exact combination of integration constants (Ci, C 2 ) that can be identified as the cosmological constant. 
The weak-curvature expansion (58) suggests that C\ ^ 0, is the crucial constant in order to attain Einstein-Hilbert 
term, but C 2 can be disposed. 

Finally we look at the wormhole's line element 

r\ 4 dr 2 



ds 2 = - j-j dt 2 + + r 2 (d6 2 + sin 2 0dcp 2 ) (59) 

and check the conditions which b(r) = ^ and $ (r) = 2 In (j^ must satisfy. The first condition which states that 



b b 2 r > is satisfied since 

b-b'r 2r „ 

— = ^>°" (60) 

The second condition implies that b' (ro) < 1. This in turn reads — r < 1 and is obvious. Finally 1 — - > means 

2 

1 — -% > and is also satisfied for r > r . The condition on $ (r) implies that it must be finite everywhere without 
root for r > ro which is trivially satisfied (except for r — > 00). To complete our example we give the embedded surface 

as z (r) — r cosh -1 . This is shown in Fig. 1 for r = 1. 



IV. CONCLUSION 



One of the most challenging problems in wormhole physics is to find an acceptable energy-momentum that will 
provide the outward push against gravitational collapse. In Einstein's general relativity, i.e. f(R) — R, this has 
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not been possible. Now, with the advent of modified theories, namely the f(R) gravity, this may turn into reality. 
Our approach to the problem of wormholes is to employ necessary existence conditions analogous to black holes. 
The existence problem of horizon in a black hole plays the similar role for throat radius in a wormhole. Expansion 
of metric functions around the throat and substitutions into Einstein's equations derive the necessary conditions. 
From these conditions, at the lowest order, we determine the metric functions that yield a traversable wormhole. In 
this process we have assumed that / (R) admits non-zero derivatives at least to third order, i.e. 4g(r ^ 0. With 
the introduction of external matter sources into f(R) we show that the weak energy condition (WEC) (at least) is 
satisfied in the construction of a wormhole in f(R) gravity. The f(R) function is explicitly determined (Eq. (52)) in 
which the Ricci scalar satisfies ^ < R < 0. Our results apply also to the case with / without a cosmological constant 
since the latter arises automatically as a combination of integration constants. In conclusion, in a generic class of 
f(R) theories satisfying the necessary conditions for existence of traversable wormhole, solutions can be supported by 
physical stress-energy tensor. The embedding diagram for the obtained wormhole is shown in Fig. 1. 
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Figure Caption: 

Fig. 1: The embedding surface (i.e. dr 2 + dz 2 = - d ^r y an d # = § ) f° r b (r) = ^ and ±z (r) = r a cosh" 1 (j^j f° r 
r = 1. 
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